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Abstract 

We calculate wide distance connected correlators in non-gaussian or- 
thogonal, unitary and symplectic random matrix ensembles by solving 
the loop equation in the 1/iV-expansion. The multi-level correlator is 
shown to be universal in large N limit. We show the algorithm to obtain 
the connected correlator to an arbitrary order in the 1/iV-expansion. 
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1 Introduction 



Recently, the universality of large iV (wide) connected correlator in random matrix 
theories is shown by Ambj0rn, Jurkiewicz and Makeenko [I]. They solved the loop 
equation for the correlators with the motivation of studying the two dimensional 
gravity. The function form of the wide correlator does not depend on the details of 
the probability distribution of a matrix ensemble. The universality classification in 
terms of this correlator is also done in uniatry ensembles by Ambj0rn and Akemann 
0. Brezin and Zee stressed this remarkably universal nature of random matrix 
theories from the view point of the physics in disordered systems ||. Although 
some other universal quantities in random matrix theories have been well-known 
already |4j], this new correlator can be calculated explicitly in various types of trace 
class ensembles in several approaches, and its universal nature can be observed 
manifestly ||. 

In this paper, we calculate the large N correlators in non-gaussian orthogonal 
and symplectic ensembles by solving the loop equation. Ambj0rn, Chekhov, Krist- 
jansen and Makeenko developed this method originally proposed in [p]] to calculate 
the correlator in principle to an arbitrary order in 1/iV-expansion in an arbitrary 
distribution of a unitary invariant ensemble ||. Here, we extend this method to 
orthogonal and symplectic ensembles, which are well-known as useful description 
schemes for the level statistics of time reversal non-integrable quantum systems with- 
out and with spin degrees of freedom, respectively [Q. Calculation scheme for wide 
connected correlator to an arbitrary order in the 1/iV-expansion have never been 
obtained in non-gaussian orthogonal and symplectic ensembles, though its leading 
order is calculated by solving an integral equation 0, its universality is shown by 
the super matrix method || and critical phenomena in an orthogonal ensemble have 
been studied by Bren and Neuberger as a model of two dimensional gravity on un- 
oriented surface ||. Here, we derive a loop equation of orthogonal and symplectic 
ensembles in the eigenvalue representation to calculate their wide connected correla- 
tor. We find a different nontrivial term from the unitary case in their loop equation. 
At the leading order in the 1/iV-expansion, the extension of the solving scheme is 
straightforward, while at the next to leading order crucial modification is necessary. 
Here we give a general scheme to solve the loop equation up to an arbitrary order 
of the 1/iV-expansion. 



2 The loop equation 



First we derive the loop equation for orthogonal(/5 = 1), unitary (/3 = 2) and 
symplectic {(3 = 4) invariant ensembles. Their eigenvalue distribution is given by 

P(Ai, • • • , X N ) = Z^AiXt, • ■ • , X N )fexp £ V(Ai)) , (1) 
where Z is normalization constant 

/oo poo ( N \ 

dX 1 ■ ■ ■ / d\ N \A(X 1 , • ■ • , A^exp -jV£ V(\i) , (2) 
-CO J — OO \ I / 
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Van der Monde determinant is denned by 

A(Ai, • • • , Aat) = II (A* -Ay), 

l<i<j<N 

and the potential is a polynomial V(x) = J2k=i ^k xk ■ The Jacobian of an infinitesi- 
mal transformation of the integration variables A^ = \ + gives 



d»\' |a(a;,...,a' a ,)|« 

This transformation in the integral (@) gives the loop equation 



.(3) 



(4) 

i/ \^ / H z_A i i=l Z — 

where (•*•) = / d N X(- ■ -)P(Ai, • • • , Aat). Here, we define one point green function 



1 A 1 



and /-level connected correlator 

N -i N -i 

gk ■ ■ ■ , *,) = e — V • • • E — L r >c 

i 1= l ^1 _ A h i != l ^ - A »! 

As discussed by Ambj0rn, Chekhov Kristjansen and Makeenko ||, the loop equation 
@ can be rewritten in terms of the one point green function and the connected two 
level corelator as the following closed form 

where the contour C encircles the singularities of G(z) . In the case of a finite N, these 
must consist of pole singularities, while in the 1/iV-expansion these make up cuts. 
Here we assume the coefficient in an arbitrary order has cuts (xi, x 2 ), ■ ■ ■ , (x 2s -i, x 2s ). 
The third term in the left hand side in eq@ is absence in the unitary case. This 
term requires a necessary modification of the 1/iV-expansion 

OO J 

G{zx, ••■,^) = E jfrGnfa* ■ - • > z i)- 

n=l iV 

Note that there is no odd degree term in the unitary ensemble. The coefficient G n (z) 
and G n (z, z) of an arbitrary order satisfies the following recursion relation 

n = 0, XG (z) = 0, (6) 
n = l, XG 1 (z)= (l-±\d z G (z), (7) 

n > 2, XG n (z) = E G n . m (z)G m (z) + ( 1 - |) 9,G n _i(z) + G n _ 2 (^, (8) 

m=l V ^/ 
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where we define the operator X acting on an arbitrary function f(z) by 

2 r dwV'(w) 



-/(to) - 2G„(z)/(z). 



j3 Jc 2m z — w' 

We shall derive large N correlators by solving this equation recursively. 



(9) 



3 Large N limit 

Here we obtain Gq(z), Gq(z,w) and Go(p,q,r) in a non-gaussian ensemble. In the 
leading order n = 0, the equation for Go can be written as 







G {z) 2 -V'{z)G {z)-R{z, V)=0, 



(10) 



where R(z, V) is defined by 



I N J 

i = l 1 

The quadratic equation (|Hj) is solved with respect to Go 

G (z) = 1 (V(z) - y/V'(z)* + 2/3R(z,V) 



If we assume s-cut solution, Go should have the form with the certain polynomial 
M{z). 

1 



(VoC) = ^|V"(«)-M(^ 



2s 



(12) 



This assumption gives the following integral expression for M 



M{z) 



dw 1 

Coo 2m W — Z 



M(w) 



dw V'(w) - pG {w) 



c x 2m ( w _ z )Jn-=i(^-x 



where the contour is a circle with the infinitely large radius. The second term 
of the integral vanishes due to the asymptotic behavior of Gq(z) ~ l/z for \z\ ~ oo 
and therefore 

dw V'(w) 



M(z) 



CoO 27T2 (nil 



w-z)JlliU(w-Xi) 



Substituting this into eq(12), one obtains 



Z Xi 



p l JCoo 2m w — z \ " w — Xi 



2m 



dw V'(w) 



(3 Jc 2m' w — z\f=\ 



2s 

n 



Z Xi 



W — X; 



(13) 
(14) 
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The large z asymptotics of Go in this expression gives 

dw w r V'{w) 



8 r 



r = 0, • • • , s. 



These equations and the arc independence of the chemical potencial 



(15) 







x 2k+l 



d\M(X) t 



■'■2k 



2s 



k 



(16) 



determine the location of the edges Xk, (k — 1, • • • , 2s) of the cuts |T(J. In the s — 1 
case, the edge values Xi, x 2 are determined only by the condition ([HJ). Here we 
calculate a wide two level correlator in the s — 1 case in the large N leading order. 
The green function is obtained by acting the loop insertion operator 



d _ _ ~ j d 
dV{w) = ~^ 1 ^+ T ^' 



(17) 



on the free energy F = -™ log Z 



, dF 1 
dV{z) z 

The /(> 2) level correlator is obtained as well 

d 



d 



Thanks to the formula 



dV{z x ) dV( Zl ) 
dV'(z) -1 



(19) 



dV(w) (z — w) 2 ' 

we can differentiate the one point green function G (z) given by eqflHD to calculate 
the two level correlator in the leading order of the 1/iV-expansion 



G (p,q) 



d 1 



dw V'(w) 



n 



p - Xi 



1 , 



+ 



dV(q) (5 Jc 2iri w — p\ fj[w — 

1 / dw —1 1 {p — x\){p — x 2 ) 
(3 Jc 2iri (w — q) 2 p — w\ (w — X\)(w — x 2 ) 

dw V'{w) 



(P - Xl)(p -X 2 ) 1 y. 



1 



dx k 



(3 Jc 2ixi p — w \ (w — x{)(w — x 2 ) 2 ^ \ w — Xk P — Xk) dV(q) 
-1 



+ 



+ 



d 



1 



(p - xi)(p - x 2 ) 



(p- q) 2 dq\p-q\\(q- x x ) (q - x 2 ) 

1 n*-(l) IP ~~ X 2 1 A/T (l) dx 2 



_ Ml (1) , 

2 dV{q)\p-x l 



+ -M, 



lp — 


Xi 


p- 


x 2 _ 



(20) 



4 



(k) 

where the moment is defined by 



M (k) _ L dw V'{w) 



the condition ( |15D 



c 2m' ( w _ xi) k \J (w — Xi){w — x 2 ) 

and M ^ ] Mf) 



In this single arc case, Af{ ff, 1 ; and ff, 2 , is obtained only by differentiating 



M (D dxi 1 



^2 (1) #t = ; 1 (21) 



(g-x 2 ) v /(g-x 1 )(g-x 2 ) 



Substituting these into eq(pO|), we obtain the following expression of the wide two 
level correlator 

n( s pq-(x l + x 2 )(p + q)/2 + x 1 x 2 1 . . 

^oyPiQ) — 1 — ~T7 — — — yj- \ zz ) 

Pip - q) 2 \l ip - xx)(p - x 2 )(q - Xi)(q - x 2 ) P\P a > 

This result agrees with that in in the unitary case (3 = 2 and that in M for an 
arbitrary beta. We can calculate the three level correlator from this expression and 
the derivative of the edge values (|21~1) 

q f p q r \— ( X2 ~ ~ x 2)(g ~ x 2 )(r - x 2 )/M 1 (1) -jp- xi)(q - x{)(r - Xi)/M 2 (1) ] 

4/3 [(p - Xi)(q - xi){r - x x ){p - x 2 ){q - x 2 ){r - x 2 )] 3 / 2 

(23) 

This result agrees with those obtained in [TJ in the unitary case. It is obvious that 
the correlator of the level more than two takes a universal form, since the two level 
correlator depends on the distribution V(A) only through the edge values 



4 Gaussian ensemble 

Here we concern ourselves with the gaussian distribution. We calculate the one point 
green function up to second order in the 1/iV-expansion. In the gaussian model, the 
one point green function and the wide two level correlator become 

1 



G (z) = - (* - ^ 2 - 2(3) (24) 

G (Z, W) = . - —, rr. (25) 

(3{z-w) 2 ^{z 2 -2(3){w 2 -2(3) (3(z - w) 2 

The loop equation (f/]) for the next order G\{z) is calculated by an reasonable as- 
sumption that the green function has no singularity outside the contour C and the 
asymptotic behavior G\{z) ~ 1/z 2 for \z\ ~ oo give 

dw w n I \ r> ( \ 
— Gi(w) = zGi(z), 

c 2m z — w 
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which yields the special solution of the loop equation (|7|) 

G ' W 4 H) (tf^T ^) ■ (26) 

The general solution is given by adding the general solution of the homogeneous 
equation with respect to Gi(z) 

Ci c 2 z 



Vz 2 ~ W yjz l -1p 

which does not satisfy the large z asymptotics. Therefore the unique solution of 
eq(^) with the correct asymptotic behavior is given by eq([26|) . These results eq(|2~4|), 
eq(^) and eq(^) agree with those obtained by the replica trick 



The loop equation (§) for G2(z) is 

XG 2 (z) = (l - I J d.G^z) + G x {zf + G (z, z), (27) 

The two level correlator with the coincident point can be defined as a finite function 
by limiting procedure 

1 



Gq(z,z) = lim G (z, w) 



which enables us to derive the solution of eq(p7|) 

n , s (2 - (3)\2z 2 + p- 2z^^W) + 2(3 2 

G2{z) ~ - WY> 2 ' ( 

This is the unique solution of eq(p7|) as derived Gi(z). 



5 Higher orders in the 1/A^-expansion 

Now, we describe a general algorithm to solve eq(^) to an arbitrary order in non- 
gaussian case. In the orthogonal (3 = 1 and symplectic (3 = 4 cases, the derivative 
term d z G n -i(z) in the right hand side of eq(0) and eq© requires the non-trivial 
modification of the algorithm in the unitary (3 = 2 case. We have to solve the 
following type of equation for an unknown function G(z) and a given function F(z) 

XG(z) = F(z), (29) 

We assume the regularity of G(z), G (z) and F(z) outside the region encircled by 
the contour C and the large z asymptotics of Go(z) ~ 1/z and G(z) ~ 1/z 2 . A 
rewriting of the second term in the operation of X and a deformation of the contour 
C z to C and give 

G (z)G(z) = I dw G °( W ) G ( W ) = _f I + I \ dw G o(w)G(w) 



c z 2ni w — z \Jc J Coo) 2ixi w 

dw G (w)G(w) 
c 2ni z — w 



(30) 
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The contribution from the contour integral of Coo vanishes due to the large z asymp- 
totics of Gq(z) and G(z). Therefore the operation of X can be written as 



^^f* 1 ""^ 1 '-^ (31) 
Jc 2iri z — w 

In this expression of the operator X, we can show a special solution of eq(^). If 
M(z)~ l has no pole outside the region encircled the contour C, we find a special 
solution 

G s (z) = F(Z) (32) 

M(z)Juti(z - Xi) 



In general, M(z) 1 has poles outside the contour C. We have to eliminate these 
unphysical singularities from the expression. Therefore, the solution becomes 

G s (z) = F{Z) - I ^ , (33) 

2M(z)Jut 1 (z-x l ) J cuc z Am ( w _ z )M{w)Jx^L 1 {z-x i ) 



dw Fjw) 

C Am l y _ o n \M(nn\.fn? a ' 



(z-w)M(w)y/UiLi(z-Xi) 



which has no pole outside the contour C. Since the second term in the right hand 
side in eq (|33|) is annihilated by the action of the operator X with the expression 
eq(BT), this G s (z) satisfies eq(^). General solution G(z) of eq(^) is given by the 
linear combination of this special solution G s (z) and the allowed zero modes of the 
operator X. In s > 2 case, Akemann shows that the allowed zero modes are 



z° z s ~ 2 



which obey the boundary condition G(z) ~ 1/z (z — > oo) [12|. In s — 1 case, the 
boundary condition excluded any zero mode to be added to G s (z). The unique 
solution of G(z) as the green function is determined by the condition eq(^) that 
the green function which is the total derivative with respect to the loop insertion 
operator ^4 , . The equation (||) for G n (z) can be solved recursively in the above 
shown formula, and a correlator with an arbitrary level can be calculated to an 
arbitrary order in 1/iV-expansion by acting the loop insersion operator eqflT7|) on 
the one point green function G n (z). 



6 Summary 

We have calculated two level and three level correlators in the large N limit in non- 
gaussian orthogonal, unitary and symplectic ensembles. We have obtained explicit 
forms of higher orders G\(z) and G^^) in the gaussian ensemble. These results are 



consistent with those obtained in other method [11]. We have shown a new general 
algorithm to solve the loop equation recursively to calculate a correlator up to an 
arbitrary order in the 1/iV-expansion for non-gaussian orthogonal and symplectic 
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ensembles. 
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